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Abstract. A common trick for designing faster quantum adiabatic algorithms is to apply 
the adiabaticity condition locally at every instant. However it is often difficult to determine 
the instantaneous gap between the lowest two eigenvalues, which is an essential ingredient 
in the adiabaticity condition. In this paper we present a simple linear algebraic technique 
for obtaining a lower bound on the instantaneous gap even in such a situation. As an 
illustration, we investigate the adiabatic unordered search of van Dam et al. |17| and 
Roland and Cerf [15] when the non-zero entries of the diagonal final Hamiltonian are 
perturbed by a polynomial (in log N, where N is the length of the unordered list) amount. 
We use our technique to derive a bound on the running time of a local adiabatic schedule 
in terms of the minimum gap between the lowest two eigenvalues. 

1 Introduction 

Adiabatic Quantum Computation (AQC) has attracted a lot of interest in recent times. 
First introduced by Farhi et al. [TT], this paradigm of computing makes use of the adi- 
abatic theorem of quantum mechanics. Informally, the adiabatic theorem says that if 
a physical system is in the ground state of an initial Hamiltonian that evolves "slowly 
enough" to a final Hamiltonian, with a non-zero gap between the ground state and the 
first excited state of the Hamiltonian at all times, then the probability that the system 
ends up in the ground state of the final Hamiltonian approaches unity as the total time 
of evolution tends to infinity. This fact is used for solving computational problems as 
follows. To begin with, the system is in the ground state of a suitable Hamiltonian. This 
initial Hamiltonian is slowly evolved to a final Hamiltonian whose ground state repre- 
sents the solution to the problem. If the running time required for a high probability 
of reaching the ground state of the final Hamiltonian is at most polynomial in the size 
of the input, we have an efficient AQC algorithm for the problem. The maximum rate 
at which the Hamiltonian can evolve at any instant without violating the adiabaticity 
condition depends inversely on square of the gap between the two instantaneous lowest 
eigenvalues. In many cases, it is difficult to estimate the gap at every instant of the 
evolution. Then, we strike a compromise by imposing a constant global delay schedule 
determined conservatively by the minimum gap over the entire interval. However, if we 
do have an estimate of this gap for the entire duration, we can apply the adiabaticity 
condition locally and speed up the rate of evolution wherever possible. 

In the early days of AQC, efforts were largely focused on adiabatic optimization 
algorithms. Given a function <f> '■ {0, 1}™ — > M, the problem is to find an x G {0, l} n that 
minimizes 6. 



van Dam et al. [T7] and Roland and Cerf [T5] demonstrated the quantum nature of 
AQC by showing a quadratic speed-up for unordered search, matching Grover's discrete 
algorithm [12] . 

However, any adiabatic quantum computer can be efficiently simulated by a stan- 
dard quantum computer |17| . Applications of AQC techniques to small sized random 
instances of several NP complete problems like EXACT-COVER [HE], CLIQUES [6] and 
SAT }8)9|11] have been explored with some success. Aharonov et al. pQ generalized this 
model by considering loca^J final Hamiltonians instead of only diagonal ones and showed 
that this generalized model can efficiently simulate any discrete quantum computation. 
Thus, in this general sense, AQC is equivalent to discrete quantum computation. 

Adiabatic quantum computation is particularly interesting because of indications 
that it is more resilient to decoherence and implementation errors than the discrete 
model. While most schemes for implementing AQC oracles involve approximating them 
by a sequence of discrete unitary gates |3|17j . robustness of potential implementations 
of the time-dependent Hamiltonians has also been studied. For example, Childs et al. [7] 
considered errors due to environmental decoherence and imperfect implementation in 
the latter approach. If the time dependent algorithm Hamiltonian is H(t), they consid- 
ered the actual Hamiltonian to be H(t) + K(t), where K(t) is an error Hamiltonian. 
In particular, K(t) can be a perturbation in the final Hamiltonian. Through numerical 
simulations, they demonstrated robustness of AQC for small instances of combinatorial 
search problems against such errors. Aberg et al. |4)5j investigated robustness to de- 
coherence in the instantaneous eigenvectors for local [5] and global adiabatic quantum 
search [4J. They showed that as long as Hamiltonian dynamics is present, asymptotic 
time complexity is preserved in both local and global cases. However, in case of pure 
decoherence, the time complexity of local search climbs to that of classical. For global 
evolution, it becomes costlier than classical: , where N is the size of the list. 

In this paper, we show how to use simple results from linear algebra to obtain bounds 
on the running time of algorithms that obey the adiabaticity condition locally. The 
technique is useful when the gap between the two lowest eigenvalues is not known at 
every instant of the evolution. As an illustration and running example, we investigate 
the behaviour of the eigenvalue spectrum when the final oracle Hamiltonian for the 
unordered search problem is perturbed in all non-zero elements by an amount at most 
polynomial in log N. 

In the unperturbed case, the gap between the lowest two eigenvalues is specified at 
every instant by a nice closed form expression |15|17j . This makes it rather easy to ap- 
ply local adiabaticity. Perturbation deprives us of this facility. We show a work-around 
by lower-bounding the gap between the two lowest eigenvalue curves with straight lines 
whose slopes are obtained using the Wielandt-Hoffman theorem from the eigenvalue per- 
turbation theory of symmetric matrices. The schedule can then be adjusted to satisfy 
the adiabaticity condition locally. The bound obtained by our method is commensu- 
rate with existing results-we obtain only a polynomial speed-up over the global algo- 
rithm |10|16|18] . Tighter bounds on eigenvalue perturbations will possibly yield Grover 
speed-up, indicating the resilience of the adiabatic algorithm to perturbations in the final 

A Hamiltonian is said to be local if it involves interactions between only a constant number of particles. 
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Hamiltonian. However, we believe that the techniques that we introduce in this paper 
can be applied in other AQC settings as well. 

The paper is arranged as follows. The next section gives a brief discussion of the adi- 
abatic quantum computing paradigm. Section 3 discusses the adiabatic search problem 
and its perturbed version. In section 4.1 we show that for the Hamiltonian in question, 
there exists a non-zero gap between the two lowest eigenvalues at all times. The proof 
largely follows that of Rao [14] . simplified for the present case. However the minimum 
gap turns out to be exponentially small, and therefore a global schedule is of limited 
value. In section 4.2 we show our method to obtain a local schedule that provides a 
polynomial speed-up. Section 5 concludes the paper. 



2 Preliminaries 



In this section we give a brief overview of the adiabatic gap theorem and its application 
to quantum computing. For details the reader is referred to the text by Messiah [13] . 
Let H(s), (0 < s < 1), be a time dependent single-parameter Hamiltonian for a system 
having an N dimensional Hilbert space. Let the eigenstates of H(s) be given by \l;s) 
and the eigenvalues by Xi(s), with Ao(-s) < Ai(s) < . . . < Ajv-i(s) for < s < 1. Suppose 
we start with the initial state of the system 1*0(0)) as |0;0) (the ground state of H(0)) 
and apply the Hamiltonian H(s), < s < 1, to evolve it to |^(1)} at s = 1. Then the 
quantum adiabatic theorem states that for a "large enough" delay, the final state of the 
system IV'(l)) will be arbitrarily close to the ground state |0; 1) of H(l). Specifically, 
|(0; 1| , 0(1))| 2 — ► 1 if the delay schedule t(s) satisfies the adiabaticity condition at every 
s: 

{) >> g{s? ' (1) 

where g(s) is the gap between the two lowest eigenvalues at s. In case g{s) is difficult 
to determine for every s, which is the case with most Hamiltonians, we impose a more 
conservative global delay schedule for the entire duration T of the evolution: 

T > max < 5 <i||^ff(g)|| 2 

Smin 

where g m in = rnino< s <i(Ai(s) — Aq(s)). However, if we do have an estimate of g(s) for 
every s, we can use a varying delay schedule that satisfies the adiabaticity condition 
locally at every instant < s < 1. Then we can reduce the running time to 



\£H(s)\\ 2 ds 



(3) 



If H(s) is a polynomial-sized linear interpolation, ||^i?(s)||2 is a polynomial-sized quan- 
tity independent of s and the running time is of the order J s= 



=0 g(s) 
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3 Perturbed Unordered Search 



Consider a list of N = 2 n elements, say bit-strings from {0, l} n . The unordered search 
problem may be stated as follows. Given a function g : {0, l} n — * {0, 1} such that 
g(u) = for a special element u and 1 for all others, find u. While any classical al- 
gorithm requires O(N) queries to g, Grover's celebrated "discrete" quantum algorithm 
accomplishes the search in 0(y/~N) queries only [12j. A similar speed-up was demon- 
strated by van Dam et al. [17] and Roland and Cerf [15] for this problem with AQC, 
which we discuss now. To begin with, note that the initial Hamiltonian H(0) should be 
independent of the solution, with the restriction that it should not be diagonal in the 
computational basis [T7J. Moreover, the ground state of the initial Hamiltonian should 
be a uniform superposition of all candidate solutions and easy to prepare. The following 
Hamiltonian satisfies the above conditions for searching an unordered list |17|15j : 

H(0)= Y, \*)(*\> W 
^e{o,i} n \o™ 

where each \z) is a basis vector in the 'Hadamard' basis given by 

and the ground state is ^= X^ze{o i} n which is easy to construct. The final Hamil- 
tonian, then, is E ze {o,i}™\« 

In this paper, we consider the case when the final Hamiltonian is perturbed in the 
non-zero entries. In other words, the final Hamiltonian is given by 

H(l)= Yl /(*)I*>R 
ze{o,i} n 

where {\z)} form the "computational basis" of the Hilbert space of the system and 
/ : {0, 1}™ — ► M po ij / ( n ) is a function that behaves as follows. f(u) = for a special 
element u E {0, l} n . For all other elements z, f(z) > 0. The problem is to minimize /, 
that is, to find u. 

Given the initial and final Hamiltonians, define the interpolating Hamiltonian as 

H(s) = (l-s)H(0) + sH(l). (5) 

We start in the ground state of H(0) and evolve to H(l) slowly enough and end up in 
its ground state. Since f(z) is bounded by a polynomial in n for all z G {0, l} n , so is 
||JjiI(s)||2. The factor deciding the running time is therefore the denominator <7^ in . 

4 The Bounds 

In what follows, we will assume f(z) ^ f(z') for z ^ z' . The motivation for this assump- 
tion is two-fold: first, it makes for a neater presentation of our technique; secondly, given 
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the random nature of noise and implementational error, it is reasonable to assume that 
no two non-zero elements of the diagonal will be perturbed by the same amount. For 
the case when there do exist z / z' such that f(z) = f(z'), the subsequent discussion 
requires only minor modification. 



4.1 Global Evolution 



By the nature of H(s), the characteristic equation is independent of the permutations 
of the diagonal elements of the final Hamiltonian. Therefore, without loss of generality, 
we will follow the convention that f{u) = < f{z\) < . . . < f(zN-i). First, we make 
sure that the gap between the two lowest eigenvalues is non-zero at all times. For that, 
we evaluate the characteristic equation of H(s), in much the same way as in lemma 1 
of [H], and PS]. 



Lemma 1. The characteristic equation of H{s) is 



c(A) = (1 



A) 



n£i 1 (i-^+^/(^)-A) = o. 



nf=i x (i - s + sf( Zk ) - a) 



l-l 

N 



^N-l (l-s+sf(z k )-X) 
l^j=l i- 3+s f( Zj )-\ 



N 



Proof. To evaluate the eigenvalue curves, we evaluate \H{s) — \I\ = 0. Subtracting the 
last column of this determinant from all other columns and using xq for 1 — s — A, x\ for 
1 — s + sf{z\) — A and so on, up to xn-i for 1 — s + sf(zN-i) — A, we have 



Xq 







xi 

: 

-XN-l -XN-1 ■ ■ ■ 



N 
N 



Xn 



-XN-l (XN-1 



(!-«) ' 
N > 



NxN 



0. 



Expanding this determinant gives the required characteristic equation. □ 

The following lemma is reminiscent of lemma 2 of |14| , and [IB] . We provide a some- 
what simpler proof for the present case. 

Lemma 2. There exists exactly one root (i.e. an eigenvalue curve) of the characteristic 
equation in the intervals (0, 1 — s), (1 — s, 1 — s + sf{z\)), . . ., (1 — s + sf(zN-2), 1 — s + 
sf(zN-i)); for < s < I. 

Proof. For any interval, let the q(A) and c u {\) denote the value of the characteristic 
polynomial at the lower and upper boundary respectively. 

We analyze the problem as three cases: 

(1) The interval (0,1 - s): 

T\ N - 1 (\ *+*f(?,Y\ l-s s^N-i nL'i 1 (ls+sf(z k )) 
lljfc=i K*--s + sj[Zk))- -jr 2_o=i i_ s+s/ ( z .) 



In this case, q(A) = (1 — s) 



T\T IlfcLi 1 (1 — s + sf(zk)). This is a positive quantity for s G (0, 1), as can easily be 
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verified. Moreover, c u (X) = —^-jf - Y\k=i s f( z k) is negative. Therefore, there exists a root 
in the interval (0, 1 — s) for < s < 1. 
(2) The interval (1 - s, 1 - s + s/(zi)): 

Clearly, q(A) for this interval is same as c u (A) of the previous case, which is negative. 

-^n^ 1 -(/(**)-/(*!)) 



But Cu(X) = -sf(zi) 



is positive. Thus, there exists at 



least one root in this interval also. 

(3) The intervals (1 - s + sf( Zi ), l-s + sf(z i+1 )), 1 < i < N - 2: 

The values of the characteristic polynomial at the boundaries are respectively 



q(A) = -sf(zi) 



c«(A) = -sf(z i+1 ) 



^ nr=\ s(f(z k ) - f{ Zl )) ntm - /(^)) 



JV-l 



and 



¥ m=i »(/(*fc) - »(/(^) - /(^+i)) 



JV-l 



Notice that in either case, every element in the first product is negative and that in 
the second is positive. This is because f(zt) < f(zi) for 1 < k < i — 1 and f{zk) > f{zi) 
for i + l<k<N — 1. Therefore, the over-all sign is decided by the number of elements 
in the first product only. But the number of such elements in q(A) and q(A) differ by 
one. Thus q(A) and c M (A) differ in sign. Therefore, there exists a root in each of these 
intervals. 

Hence, there is a root in each interval. Given that (i) there are N open intervals 
bounded by N straight line eigenvalue curves, (ii) there is at least one root in each 
interval and (iii) there are TV roots in all, the lemma follows. □ 
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By virtue of the above lemma, we can speak of one curve being "above" another. 
We label the curves as Aq(s), Ai(s), . . . , Ajv-i(s), starting from below. Thus, Aq(s) lies 
between the lines A(s) = and Ai(s). See figure 1 for example. 

As a consequence of the above lemma, we are guaranteed a non-zero gap between 
Ao(s) and Ai(s) for all s G [0, 1]. 

What about the minimum gap between Ao(s) and Ai(s)? It turns out that this gap is 
exponentially small [18]. Making use of the properties of the present problem, we arrive 
at this conclusion in a different manner. 

We first discuss the case when f{z\) is large, say greater than 1. Consider the charac- 



ll/e=l {L-S+SJ {Zk)-A) 7TZ^j=l l-s+sf(zj)-X 



teristic polynomial c(A) = (1— s— A) 

^nti'a --+«/(**) -a). 

By lemma 2, the line A(s) = 1 — s separates Ao and Ai. Thus, if 

(i) there exists an sq where the above polynomial is divisible by two "(1 — sq — A) 
factors" and 

(ii) all other curves (A2(s) through Ajv-i(s)) are at least an inverse polynomial dis- 
tance above the point (sq, 1 — sq), 

it would imply that Ao(so) and Ai(so) are exponentially close to the line A(s) = 1 — s, 
and in turn, to each other. 

Clearly, points very close to s = are not candidates, as there are too many curves 
in the vicinity. Points close to s = 1 are also ruled out, as there is only one (1 — s — A) 
factor, corresponding to the curve Ao(s): Ai(s) is f(z\) units above. 

Consider a candidate point sq = where p(n) is a polynomial in n such that 

p(n) > maxk(f(zk))- At sq = A2 > 1 + 1 ■ Thus, by lemma 2, all curves above 
A2(s) are at least an inverse polynomial distance above 1 — so = 1 — 

Notice that the polynomial can be divided by one (1 — s — A) factor at an s such that 
rifcl^i 1 (1 — s + s f( z k) — A) — ► 0. Let us test the candidate point sq. Substituting 1 — s 
for A and yj^y for s, we get ^y^r Ilfc^i 1 ( p(n) )' wn i cn indeed tends to zero. Therefore, 
the existence of one (1 — s — A) factor at so nas been established. Let us now see if 
there exists another. On dividing the polynomial by (1 — sq — A) at sq = ^ns, we obtain 



c(A) 



tJV-1 f{z k ) 



nr 

llfc=l "^J" Np{n) 2sj=l JTJD 



, which tends to zero. 



p(n) 

Suppose now that f(zi) is "small". If f(z\) is a constant independent of n or even 
0( poiy(n) )' ^ ne argument given above holds exactly. However, if f(z\) is greater than zero 
by only an exponentially small quantity, the two (1 — s — A) factors that we obtained in 
the above discussion could correspond to A2(s), thus leaving the possibility of Ao(s) being 
at a larger distance below A(s) = 1 — s. But this is not the case. Since by assumption 
A2(s) is exponentially close to the line A(s) = 1 — s, it can be factored out from the 
polynomial, leaving two factors that correspond to Ao(s) and Ai(s). 

Thus, it turns out that the minimum gap is too small to yield any significant speed-up 
for a constant rate global adiabatic algorithm. 
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4.2 Local Evolution 



We will now obtain an upper bound on J g=0 where g(s) = Ai(s) — Ao(s), for lo- 
cal evolution. To that end, the following result from the perturbation theory of real 
symmetric matrices will be useful. 

Theorem 1. Wielandt- Hoffman Theorem (WHT): Let A and E be real symmetric ma- 
trices and let A = A + E. Let the eigenvalues of A be Xq < X\ < . . . < \n-i and those 
of A be Ao < Ai < . . . < Atv_i. 
Then 

N-l N-1N-1 
j=0 i=0 j=Q 

□ 

Let E of WHT be the 'perturbation matrix' H(s + ds) — H(s). This provides us with 
the r.h.s. of (|6|): 

k=l 

For the l.h.s., note that X\(s) . . . Av-i(-s) are N — 1 non-crossing curves packed inside a 
polynomially bounded gap and each curve is bounded by two straight lines by lemma 2. 

Thus, there are at most a polynomial q(n) number of the gaps f(zi) — that 
are at least inverse polynomially wide. All other curves are packed between the straight 
lines A = (f(zi) -l)s + l and A = - l)s + 1 where /(*) - = 0(f ) for 

a constant c. Therefore the slopes of these curves can be approximated by the slopes of 
one of the enclosing lines. Hence, the l.h.s. is 

d\i + d\\ + . . . + d\\ {n) + y, - 1 ) 2 ( ds ) 2 

fce{i,...,AT-i}\{ii,...,i g{n) } 

Substituting in El we get an upper bound on dX?, + dA? + . . . + dXf . In particular, this 

U £1 t q(n) 

is also an upper bound on ^ and This can be used to obtain an upper bound on 
the total time required, as we will see soon. 

We illustrate the technique for the case when q{n) = 1 and f(z\) > 1. 

Lemma 3. If q(n) = 1 and f(zi) > 1, 

2 A similar bound can be obtained using a general result of Ambainis and Regev ([2], Lemma 4.1). 
Nevertheless, we presented a different way to demonstrate the possibility of problem specific approaches 
which can yield tighter bounds than the general one. 
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Proof. By the preceding argument, the l.h.s. ofE]is given by 

Ylk=i(f( z k) ~ l) 2 {ds) 2 + (dX (s)) 2 + (tiAi(s)) 2 . Substituting in equation^! we get 
Ek=i(f(^) ~ 1) W + (dXo(s)f + (dX 1 (s)) 2 

< (^) W + ^(ds) 2 + - ^) W- 

Rearranging some terms we get, 

(f 1 ) 2 + (^) 2 < (TV 1 ) 2 + V + M*N-l) ~ 

+ Ek=i (if - ^v 1 ) 2 - (/(*)« ~ V) 2 ) 2 - 

Or, 

fd\ Q \2 i ^rfAi\2 ^ ( N-\ \2 i jV-1 i (f( v \ B i 1 r^ 2 ^ x •)//, 1 2N-l \ 

Simplifying and ignoring small terms, we get 



N— 2 



fc=i 



□ 



The lemma implies that |^|, |^| < ^ (/(*;v-i) - l) 2 + 2n + | ££~i 2 /S)- 

Using this, we will estimate g(s) = X\(s) — Xq(s). We conservatively approximate 
Xq(s) and Ai(s) by straight lines to get a lower bound on g(s) = X\(s) — Xo(s). 

We divide the interval [0,1] into three parts [0, a], [a, b] and [6,1], corresponding to 
the parts when Ao approaches A = 1 — s, when both Ao and Ai are close to A = 1 — s, 
and when Ai(s) rises away from A = 1 — s respectively. 

Consider the first interval. Suppose Iq(s) and l\(s) are lines such that Iq(s) > Xq(s) 
and l\(s) < Ai(s). Then the gap between these lines is a lower bound on Ai(s) — Ao(s) 
in the interval (0, a). 

Denote \J (/(.z/v-i) — l) 2 + 2n + jjj X]£i 2 /( z fc) — 1 by m. We take = ms and 
/J(s) = 1 — s. Thus, the gap in the first interval is given by g%(s) > 1 — s — ms. Similarly, 
for the third interval, we choose Iq(s) = 1 — s and if (s) = ms — m + f(zi). Therefore, the 
gap 53(5) in this interval is greater than ms — m + f(z\) — (1 — s). We will conservatively 
take the gap in the entire second interval to be g m in- See figure 2 for a rough sketch. 

Let a be specified by a point on Iq(s) = ms that lies vertically below the line A = 1 — s 
by a distance g m % n - Thus, a = -^ff 1 ■ Similarly, b is specified by a point on the line 
lf( s ) = ms -m + f{zi) that is above A = 1 - s by g min . Then, 6 = m+1 ~^ 1 1 )+gmin . 

For the gap g%(s), recall that one of Ao and Ai is away from the line 1 — s by a margin 
of drain- Since this gives the deviation of only one of Ao and Ai from 1 — s, this is a 
conservative estimate of 52 (s)- 

Therefore the total delay factor is given by 

ds f b ds f 1 ds 



9i (s) 2 J a 92{s) 2 J b g 3 (s) 2 

f a ds 1 f b , f 1 ds 
< / rj + ~2 — ds+ 

J ° (l-s(l+m)J 9minJa Jb + ]_)( s — 1) + f (z\) 
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Fig. 2. Bold lines show actual eigenvalue curves and the dashed lines show lines that are 
assumed to bound the gap. 

_ 1 /_L j\ m + 2g min - f(zi) 1 /_L 1 \ 

m + 1 \g min J 9 m in( m + 1 ) m + l \g m in /Oi) ' 
m - f(zi) 

~ (™ + 1 )i„' 

Let us summarize: 

Theorem 2. Let the non-zero elements of the final Hamiltonian of the local adiabatic 
search algorithm be f(z\) / f{z2) / f(zN-i), where each f(zi) is of size at most 
O(logiV). Then, the time taken to evolve to the solution state of the final Hamilto- 
nian is 0{D where m = yj (/(zjv-i) - l) 2 + 2n + ^ T,k=i f( z k) ~ 1, 9min = 
min < s <i(Ai(s) - X (s)), and D = \\£H(s)\\ 2 . 

To obtain tighter upper bounds on the running time, we would require the interval 
b — a to be shorter: short enough to balance the denominator o 2 • in the second term. 
This in turn requires tighter bounds on and 

5 Conclusions 

We introduced a technique for obtaining upper bounds on the running time of adiabatic 
quantum algorithms if the eigenvalue spectrum behaves in certain ways. We used this 
technique to investigate the robustness of the adiabatic quantum algorithm for unordered 
search when the final Hamiltonian is perturbed in the non-zero entries. Interesting open 
problems include tightening of the bounds, and application of our technique to other 
quantum adiabatic algorithms. 
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